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We study various improved staggered quark Dirac operators on quenched gluon backgrounds in 
lattice QCD generated using a Symanzik-improved gluon action. We find a clear separation of the 
spectrum into would-be zero modes and others. The number of would-be zero modes depends on 
the topological charge as expected from the Index Theorem, and their chirality expectation value 
is large (« 0.7). The remaining modes have low chirality and show clear signs of clustering into 
quartets and approaching the random matrix theory predictions for all topological charge sectors. 
We conclude that improvement of the fermionic and gauge actions moves the staggered quarks closer 
to the continuum limit where they respond correctly to QCD topology. 
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INTRODUCTION 

It has been widely held that lattice staggered quarks 
are insensitive to the topology of the underlying gauge 
fields. The low-lying spectrum of the Dirac operator has 
neither shown the number of chiral (near-) zero modes 
anticipated from the Index Theorem, nor have the eigen- 
values lain on the expected universal distributions. 

Here we present evidence that this is not a generic 
failing of staggered quarks, but simply a problem of dis- 
cretisation errors, and that the use of improved staggered 
Dirac operators clarifies the situation and points to the 
correct continuum behaviour. This allows the use of im- 
proved staggered quarks in lattice QCD to study topo- 
logically sensitive states, such as those associated with 
the axial anomaly (principally the rj' meson). It also has 
a bearing on establishing the effect of taking the fourth 
root of the staggered determinant to represent one flavour 
of staggered sea quarks. 

We begin by reviewing our understanding in the con- 
tinuum. The eigenmodes of the (anti-Hermitian, gauge 
covariant) massless Dirac operator are given by 



pfs = iXsfs 



Xse 



(1) 



where we use orthonormalised eigenvectors, fjft ~ Sst- 
As {^,75} = 0, the spectrum is symmetric about zero: 
if Xs ^ 0, then 75/s is also an eigenvector with eigenvalue 
— iAs, and chirality Xs = /1 75/s = 0- The zero modes. 
As = 0, can be chosen with definite chirality: Xs = il- In 
general there are n± such modes, whose relative number 
is fixed by the (gluonic) topological charge 
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via the Atiyah-Singer Index Theorem [l|, |2| : 



Q — mTi 



75 



n+ 



(2) 



(3) 



where m is the quark mass. Based on 3] , it has been sug- 
gested that (for a sufhciently large volume) the non-zero 
low-lying eigenmodes take values from a universal distri- 
bution y scaled by a QCD-specific quantity (the chiral 
condensate). The universality class is determined by the 
chiral symmetries of QCD with separate predictions for 
each sector of fixed topological charge. The distributions 
can be derived from any theory in the correct universal- 
ity class, such as ensembles of random matrices [5|,|a| (for 
a review of other theories, see |7|). 

The four-dimensional lattice staggered massless Dirac 
operator is anti-Hermitian, with a purely imaginary spec- 
trum i\s. It represents Nt — A "tastes" of fermions that 
interact via highly-virtual gluon exchange at finite lattice 
spacing, causing taste-symmetry violations [g]. These 
vanish in the continuum limit (as a?) and we then expect 
to recover a four-fold degeneracy in the spectrum. A 
remnant of continuum chiral symmetry in the staggered 
action gives a local and taste-non-singlet 75 operator that 
guarantees that the spectrum is symmetric about zero, 
as in the continuum. The 75 operator, 7^'', relevant to 
the index theorem must be a taste-singlet one, however, 
since only this can couple to the vacuum correctly |9(. 
As {p, 75^*} ^ 0, there is no exact Index Theorem and 
all eigenmodes in principle contribute to Eqn. ^ 10] . 

If the gauge field is sufficiently close to the continuum 
limit, however, we expect to see the continuum features 
developing. There should be 2|Q| near-zero modes on ei- 
ther side of zero, whose chiralities are close to unity. The 
taste-singlet 75'' operator is not conserved so we expect 
a renormalisation to achieve a value of 1 in the contin- 
uum. The remaining modes should have chirality near 
zero, and come in approximately degenerate quartets on 
either side of zero. The values of the eigenvalue quartets 
should be described by the same universal distribution 
as continuum QCD, up to a renormalisation of the chiral 
condensate. 



For thermalised lattices at finite lattice spacing, the 
fluctuations in the gauge field can lead to a breakdown 
of this picture. This happens for the simplest, "one-link" 
(naive) staggered operator with the unimproved, Wil- 
son gauge action at values of the gauge coupling used 
in present-day simulations. The breakdown is seen in 
two ways. Firstly, there is no clear separation (in eigen- 
value or chirality) of near-zero modes of topological origin 
from the remainder of the spectrum [ifll, |ll|, |l3, llj, llJ| ■ 
In addition, the low-lying eigenvalues do not follow the 
predictions from universality. In fact, the eigenvalues in 
each sector of charge Q ^ all follow the distribution cor- 
responding to the sector with = [TsUllllTllli Hgj . 
(It should be noted that the method we follow here, of 
grouping the eigenvalues into quartets, was not followed 
because this feature of the spectrum was not evident.) 

This failure can be ascribed to taste-changing interac- 
tions and the lack of a good continuum chiral symmetry 
in the one-link staggered Dirac operator: good agree- 
ment with predictions for all topological charge sectors 
has been seen for Dirac operators obey ing the Ginsparg- 
Wilson relation [13, El Hi HI H if . 

Over the past few years considerable advances have 
been made in lattice QCD phenomenology through the 
use of so-called "improved staggered" fermion formula- 
tions [26|. The goal of this programme is to systemati- 
cally reduce the lattice artefact taste interactions, and 
it is reasonable to expect that this will improve the 
continuum-like chiral properties of the fermions. The 
variation of the topological susceptibility with the sea 
quark mass, in particular, shows that this is so: whilst it 
was largely insensitive to the presence of one-link stag- 
gered sea quarks J23, 123, 12^, U^ , the improved operator 
gives a variation with m that agrees well with theoret- 
ical expectations [Sfl |31| • It is thus pertinent to study 
in more detail the spectrum of the improved staggered 
Dirac operator, and we shall show here that this also 
now gives signs of converging to the expected results, in 
contrast to previous studies. 



RESULTS 

On a Euclidean lattice with lattice spacing a, the one- 
link massless staggered operator is 



I 4 



(4) 



Al=l 



with i]i = I and rj^ = (— l)'^''7yi,_i;_ We study also three 
improved operators: the ASQTAD [J. Is^ Isj. IJ^I . the hyp 
I and the Highly Improved Staggered Quark (hisq) 
ll, ISol . These operators use "smeared" gauge fields in 
place of the U field above, obtained by multiplying U 
fields along combinations of bent paths from the start 
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FIG. 1: The positive half of a typical low- lying eigenmode 
spectrum for a configuration of Q = 2, for various staggered 
fermion formulations. The bottom panel gives the absolute 
value of the eigenvalue. As, ordered according to increasing 
size. The x axis is then simply eigenvalue number. The top 
panel is the chirality of the modes. The hyp action gives 
results very similar to HiSQ and is not plotted for clarity. 



to end points of the original link. This reduces the cou- 
pling to highly virtual gluons and suppresses the taste- 
changing interactions. The ASQTAD action uses a "FAT7" 
smearing, which includes paths made of up to seven links, 
and tadpole improvement. The hyp action uses an hyper- 
cubic blocking procedure, involving rcunitarization back 
onto SU(3), and the HiSQ action uses two applications 
of the FAT7 smearing, and also includes rcunitarization. 
The taste-changing interactions are most suppressed for 
the HYP and HISQ cases. The HISQ and ASQTAD cases are 
completely a^ improved at leading order by additional 
terms which correct for errors in the simple derivative 
above. 

We calculate the eigenmodes for these Dirac operators 
on an ensemble of quenched (no sea quarks) SU(3) gauge 
configurations. The gauge action is Symanzik-improved 
at tree-level with tadpole improvement so that remain- 
ing discretisation errors from the gluon field are a small 
number times ajO^. The lattice spacing is 0.093 fm, rep- 
resenting a standard ensemble in present-day lattice sim- 
ulations [37J. The majority of our results are from 1000 
configs on a periodic lattice with 16"* sites, which should 
be large enough that finite volume effects on the low- 
lying spectrum are negligible |25|. We have also studied 
a larger volume, 24^. On each configuration we deter- 
mine the topological charge by two standard methods 
that involve cooling the gauge fields [31| . The calculation 
of gluonic topological charge always involves ambiguities 
and we discard those configurations (10% of the ensem- 
ble) for which the two cooling methods do not agree, to 
leave a sample for which we are confident that we have a 
robust estimate of Q. We stress at this point that coohng 
is used solely to determine the gluonic topological charge, 
and that all the Dirac eigenmode calculations are carried 
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FIG. 2: A scatter plot for difFerent staggered quark formu- 
lations, with absolute value of the chirality on the y axis and 
eigenvalue As on the x axis. The lowest 50 eigenvalues for 147 
configurations are plotted. 



out on fully thcrmalised, "raw" configurations. 

Fig. n compares the low-lying modes of the spectra 
of the various staggered quark formulations on a typical 
background with Q = 2. We show the eigenvalue As and 
chirality for the upper half of the spectrum. We see that 
the Index Theorem is well approximated for the more 
improved quark formulation. Specifically, there is a clear 
delineation between the near-zero modes (small eigen- 
values and large chirality) and the rest of the spectrum. 
The number of near-zero modes is 2|(5|, as expected. The 
other eigenvalues have small chirality, and divide clearly 
into quartets. The taste-singlet 75* operator used for the 
chirality is a point-split 4-link operator. We measure it 
here by inserting U fields to make it gauge-invariant. We 
will report elsewhere on the dependence of the chirality 
on the operator used. 

As the near-zero modes for the more improved actions 
clearly separate and have well defined chirality, we may 
define an index Q — n_|_ — n_ on each configuration. Q is 
then strongly correlated with Q. If we count eigenmodes 
with absolute value of the chirality above 0.65 in n±, for 
example, we find that in about 90 % of configurations 
Q and Q are the same, satisfying the Index Theorem. 
Indeed, measuring Q from the Index Theorem then be- 
comes as reliable as measuring it from gluonic methods. 

In Fig. 13 we show a scatter plot of the absolute value 
of the chirality versus the (absolute value of) eigenvalues, 
for different operators. We can see the formation of a gap 
between modes of small and large chirality as we improve 
the staggered operators, as well as the overall increase in 



chirality of the large chirality modes. 

We turn now to the non-zero modes. Subtracting the 
2\Q\ near-zero modes from the spectrum, we group the 
other eigenvalues, ordered by size, into sets of four, as in- 
dicated in the HISQ and hyp cases in Fig. ^ We call the 
average of these sets Ai.2,.... In Fig. Owe plot the ratios 
(As)Q/(At)Q (denoted by '^s/t"), where the expectation 
values (•)q are over configurations with gluonic topolog- 
ical charges ±Q only. Also shown are the universal pre- 
dictions. There is a clear dependence of the ratios on Q, 
in marked contrast with previous results, which showed 
a precise agreement with the Q = predictions for all 
sectors. 

The results are systematically slightly lower than the 
theoretical predictions, especially for the ratios involv- 
ing higher eigenvalues. This would be consistent with 
finite volume effects as in |25j. There is also a small 
but systematic difference between the one-link and the 
improved actions, with the improved results showing a 
better agreement with the theoretical values. 

An important point to make here is that it is neces- 
sary to group the eigenvalues as explained above to get 
sensible results. If one ignores the near zero modes, or 
does not group in quartets, ratios which are close to one 
or very large will result. This is strong evidence that the 
four tastes are showing up in the spectrum, even where 
it is not directly evident in the spectrum itself. We are 
undertaking further analysis to understand why previ- 
ous results with naive staggered quarks on unimproved 
giuon fields were in agreement with universal distribu- 
tions for Q — for all Q values when this procedure was 
not followed. It seems clear from our analysis here that 
results on unimproved gluon fields with a fine enough 
lattice spacing should show Q-dependent universal dis- 
tributions when zero modes and quartets are taken into 
account. Preliminary results confirm this expectation. 



CONCLUSIONS AND OUTLOOK 

Improved staggered fermions are not blind to the topol- 
ogy, but in fact reproduce well the predictions of the In- 
dex Theorem, and the universality of ratios of eigenvalues 
as a function of topological sector. This means we can 
have confidence in using them to attack the questions 
arising from the axial anomaly in QCD. 

We also remark that the fact that the 4-fold taste de- 
generacy of staggered quarks is becoming clear in the 
spectrum is encouraging for the programme of establish- 
ing the eflFect of taking the fourth root of the staggered de- 
terminant to represent one flavor of staggered sea quarks. 
This programme requires an analysis in the taste basis 
and progress towards this is now possible. 

More extensive studies of finite volume and lattice 
spacing effects and analysis of the eigenvectors are un- 
derway and will be reported elsewhere. 
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FIG. 3: The ratios of expectation values of small eigenvalues 
(see text for notation) compared with the predictions based 
on a universal distribution (horizontal lines) for topological 
charge sectors 0, 1 and 2. 



In the later stages of this study we became aware of 
work on a related topic ^38] . 
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